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Abstract 

Wo  ooiisUlor  llio  form  ol  i.  Iivoakiui!  wavo  <li.s.si|>aUoi.  toriii  for  iiac  “ 

stodiaslir  wave  evolutiou  iiioaols.  A  limio.loiuai.i  lloiissiiioscl  inmlol  is  “f 

;us  suf'^esUHl  bv  EMeberky  and  Hatties  (1990). 

Introduction 

uooo,a,,  11,0.  ,.avo  I 

wave  dissipation  term  m  spcctriil  .  ‘  .  (Krii.nt  -us  a  function  of  fre- 
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related  to  surface  displacement  i;  accoidin^,  tc 
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where  tlie  omitted  terms  are  related  to  shoaling  and  nonlinear  interaction  effects. 
Masc  and  Kirby  (1992)  i)rt>|)o.sed  a  lorm  lor  <v„  given  by 


with 


—  <»()  +  (/n//)^Ol 


(3) 


rvo  =  Fft 


n,  =:(1-F)/^ 


(4) 

(3) 


and  wh(M  (!  (i  is  deUn  inined  from  a  bulk  dissipation  model  such  ;is  the  one  of  Thornton 
and  Guza  (1983).  Based  on  analysis  of  laboratory  data  and  numerical  results,  Mase 
and  Kilby  (.hose  to  set  F  —  0.5,  indicating  a  dissipation  term  with  a  partial  dependence 
on  the  s(|uaie  of  the  frocpiency.  They  also  found  that  choosing  F  =  0.0,  corresponding 
loan  /'^  d(>pemlem:e  for  tiie  entire  di.ssipation  term,  destroyed  the  tail  of  the  computed 
power  spc'ctrum  in  very  sliallow  water  but  had  little  impact  on  the  evolution  of  spectral 
shape  away  from  the  shallowest  measuring  gages. 

Eldeberky  ami  Battjes  (1990)  have  suggested  that  a  similar  formulation  corre¬ 
sponding  to  tli<‘  (  lioire  F  ~  1.0,  spreading  the  «lis.si])ation  term  uniformly  over  all 
frcquenc.ie.s,  should  lx;  utilized,  and  slnnved  that  an  aderpiate  de.scriirtion  of  power 
spectrum  evolution  was  obtained  in  .several  simulations  of  held  data.  Eldeberky  and 
Battjes  did  not  consiiler  the  efiect  ol  this  choice  on  the  (>volution  of  higher  statistical 
moincnt.s.  Moie  tecmitly.  C/hen  c't  al  (1990j  referiasl  to  a.s  GGE)  havi?  examined  a 
numbrn-  ol  laboratory  and  fiehl  cases.  Th(>y  have  shown  that  the  (estimates  of  jrower 
spectrum  (Solution  are  relatively  in.sensitive  to  the  choice  of  F,  with  error-minimizing 
F  value, s  (jccupying  tin;  entire  range  0  <  F  <  1  for  various  firdd  and  laboratory  cases. 
Aside  from  tin?  Ma.s(.*  and  Kirby  cjuse,  error  measures  changed  by  as  little  as  20%  over 
the  entire  rangi;  of  vabu's.  In  contrttst,  all  data  .sots  support  the  choice  of  F  =  0.0  when 
****^‘*^*' I e.s  b<i.s(*d  on  thir<l-moment  statistics  are  introduced,  with  the  exception 
of  the  Ma.se-Kirby  data  set.  (This  last  discrepancy  is  fairly  weak,  however).  There  is 
a  clear  trend  towards  increasing  error  with  increasing  F  in  most  datasets. 

In  this  talk,  the  problem  of  determining  the  form  of  the  spectral  dissipation  term  is 
approac.hed  fiom  a  different  direction.  Instead  of  considering  a  bulk  energy  decay  and 
an  arbitrary  distribution  of  di.s.sipation  over  /,  we  instead  consider  the  structure  of 
the  dissipation  term  in  a  time-d(*pendent  Doussine.s(]  model  setting,  and  consider  the 
contribution  of  the  term  to  energy  lo.s.s  and  the  structure  of  that  loss  in  the  frequency 
domiiin.  Ihis  lo.s.s  is  then  rtdated  to  the  spectral  evolution  equations,  and  the  form  of 
the  dissipation  t(U'm  is  deduccsl.  We  find  that,  in  general,  the  dissipation  coefficient 
is  distributed  moK*  or  le.ss  a.s  \  wln.ue  5,,(/)  is  the  power  spectrum  of  the 

surface  di.splacement  j/.  for  the  r  a.se  of  a  smooth  spectrum,  this  residt  indicates  an 
/*  dependence  in  th(‘  dis.sipation  coefficient  in  the  surfzone,  consistent  with  CGE’s 
results  with  /•  --  ().().  1  he  (  I inclusions  h(!i(r  are  ba.scsl  on  an  examination  of  the  Mas(! 
and  Kilby  ( 1992)  data  set,  which  is  well  imxh'lh'd  by  the  chosen  time-domain  breaking 
wave  model. 
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Tinie-Doniain  Model  lor  Wave  Breaking 
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I'lie  lileratuic  includes  examples  of  a  number  of  lormulations  for  computing  wiv^ 
brenkiiig  in  the  ront(‘Xt  cf  BoiissinesM  wave  models.  Three  examples  include  an  eddy* 
vuscosify  model  (Zell,  191)1),  asurlace  roller  model  with  an  iissumcd  two-layer  veloc^W 
profile  (Schaffer  et  al,  1993),  and  a  surface  roller  model  with  a  computed  honzontalj 
profile  (Svends(‘i.  <‘t  al,  1990).  Schematically,  the  eddy  vi.scosity  model  may  be  written- 

as 

II,  -p  un^  +  IJ1U  +  dispersive  terms  -  =  0 

The  simple  roller  model  of  Schaffer  et  al  (1993)  is  written  in  the  context  of  a  mod^ 
for  total  volume  flux,  and  may  be  written  schematically  as  ^ 

p  ( P'^/lA  +  qlltlx  +  dispersive  terms  +  /Zr  =  0  (7^ 

where 

/?  =  /  ^^  («?ol  -  W?rrol) 

4 

The  more  <-omplex  rolhrr  model  of  Sveml.sen  et  al  (1990)  is  essentially  of  the  8ame| 
form.  While  these  models  ilifler  in  both  form  and  theoretical  intent,  it  may  be  showa^ 
that  the  numerical  rei)resentations  of  each  <li.s.sii)ative  term  are  similar.  In  particular? 
the  contribution  of  each  dissipation  term  is  highly  localized  in  space  and  time,  since,; 
it  is  concentrated  on  the  front  face  of  the  breaking  wave.  Further,  the  contribution^Ofj 
each  dis.sii)ation  term  is  about  the  same  size,  since  each  successfully  calibrated  moddi 
must  extract  the  same  amount  of  energy.  An  illustration  <»f  this  fact  is  given  by  Figurt 

9  of  Sveiulsmi  (*t  al  (1990).  r  i  ’^1 

Sinc<*  our  primary  goal  1m*Iow  is  to  examine  the  spirtral  signature  of  the  wavjl 

energy  decay,  we  can  conclmle  that,  due  to  the  structural  similarity  of  the  vario^ 
breaking  moilels  in  the  time  domain,  it  should  not  matter  which  of  the  existing  mod-j 
els  is  used  to  perform  the  analysis.  The  analysis  will  be  based  on  the  eddy  viscosit^^ 
model  of  Zelt  (1991)  for  two  reasons:  it  is  already  incorporated  in  an  existing  tim^ 
domain  Doussinesi,  model  (Wei  and  Kirby,  1990),  and  it  is  .simple  to  interpret  the 
terms  in  the  eddy  vi.scosity  ino.lel  in  terms  of  measured  sea  surface  elevations,  as  d^ 

scribed  l)elow. 


//  =  /*  +  V 


Lentling-Ortlcr  Energy  Unlnnce 

,.  =  it/jj  oi  repre.seiit  tlie  breaking-imluce*!  momentum  <leficit  per  dm| 

depth.  Then,  eai  li  of  the  models  above  may  be  written  in  the  form  ^ 

n,  +  uiir  -f  (Jijr  +  ilispersive  terms  -1  iv  =  d  (10)| 

Neglecting  noidiuear  and  dispersive  effects,  we  have 

ui  d-  !JVx  +  =  d 

=  (»  (12)| 


0 

0 


tiplying  (11)  l>y  (12)  by  p<jv  ^“<1  adding  then  gives 


^  E  =  ^-pgri"  + 

t  the  local  energy  density/unit  surface  area, 

F  =  pghui) 

b  the  flux  of  energy  in  the  x  direction,  and  ^ 

fj,  =  phur^ 

i  the  local  rate  of  energy  decay.  Each  of  these  quantities  may  be  averaged  in  time, 
■yielding  (for  a  stationary  wave  process) 

{F)r  =  -(ffc) 

The  average  or  bulk  energy  decay  (a)  can  be  specified  according  to  models  such  as  the 
of  Thornton  and  Gu/a  (1983),  which  models  the  data  set  considered  below  qm 
well.  Each  of  the  quantities  in  (17)  may  be  thought  of  as  the  sum  of  contributions 
from  ciich  frc(iu(Micy  to  tlic  total  value,  i.c,, 


I 


(F)  =  ^F,.;  (f6)  =  X^ffcr« 


We  will  attach  a  meaning  to  each  of  these  component  terms  in 

In  order  to  proceed  further,  we  need  to  choose  a  model  to  evaluate  ((..  This 

be  done  using  the  Zelt  (1991)  eddy  viscosity  model. 

Eddy  Viscosity  Model 

The  eddy  viscosity  appearing  in  (6)  is  written  by  Zelt  (1991)  as 

Uf,  =  i  =  By  {h  +  Tj)  (13) 

hpre  -Y  =  2  is  a  mixing  length  parameter  determined  by  Heitner  Housner  (1970) 
ind  chlscn  so  that  the  resulting  model  correctly  predicts  the  width  of  n  hydraulic 

jump.  The  factor  D  is  given  by 


B  =  {  U:r/K  -  1 
n 


Ujc  <  2u* 

2ii’  <  iix  <  K 

Ur  >  llx 


„d  provides  a  somewhat  smooth, si  on.set  of  breaking  dissipation  when  I''.-- 

itg  criterion  is  excecled.  The  breaking  criterion  is  given  m  terms  of  a  critical  velocity 

divergence,  chosen  to  be  - 


traLrbcl«7  "•  waves  will-  be.£ 

cd  btlou  Tl.e  inodel  .s  particularly  useful  to  us  in  the  present  context  becaii 

the  terms  in  the  model  may  be  evaluated  (to  leading  order)  using  time  derivative 
the  surface  elevation.  Using  ^  -j/./U,  we  get  r'vaiives 

n  =  Dy'^h^Uj,  ss  Dy'^hi)i 

\  1  :  >  27;; 

^  =  j  Vi/v:  -  1  ;  7/;  <  in  <  27;; 

I  ;  m  <  Vt 

7;,*  =  0.3  v/^ 

Finally,  the  instantaneous  energy  dissipation  may  be  written  as 


ffc  =  -phu  ~  -P 


The  advantage  of  this  formulation  is  clear  in  the  context  of  evaluating  expenmen 
results,  since  the  energy  lo.ss  term  that  would  be  predicted  by  the  numerical  mo3 

tion  effects,  W^e  may  proceed  without  actually  running  the  model  in  the  majoriJ 
es,  provided  that  the  model  as  formulated  i.s  known  to  be  an  accurate  predictof 
the  wave  field  in  sample  representative  case.s. 

Laboratory  Data 

119Q^*r'  taken  from  Run  2  of  Mase  and  Kirbl 

{  -).  le  present  results  have  ln>en  reproduced  for  a  number  of  other  data  set^ 

and  a  more  comprehensive  view  of  the  study  will  l,e  publisheil  elsewhere.  The  sinjS 

sttdiedTdtt^^^  ""  ^  of  conditi^l 

Figure  1  shows  a  schemalic  of  the  experimental  facility.  The  experimental  waw 

ZtT"  ''I  "  ^‘‘^--‘•^^‘^-kowitz  spectrum  generate<l  in  47cm  of  wat^ 

with  a  i)cak  frequency  /  =  l//.  and  a  significant  wave  height  of  Gem.  Waves 

D-i^'for  tr*"^  capacitance  wave  gages  at  twelve  stations  across  the  1:20  beach  profU^I 
ata  hr  the  analysis  below  i.s  taken  from  the  iiieiisurement  at  the  h  =  10cm  deptlf 
This  depth  correspomls  t<.  a  point  where  the  probability  of  breaking  is  increasing 
apidly  but  the  saturate<l  inner  surf/one  lues  not  yet  been  established.  -t; 

at  /,  -"in'"  nieiLsured  and  comput(>d  time  series  of  surface  elevatiof 

at  h  10cm  IS  shown  m  Figure  2,  which  indicates  an  accurate  reproduction  of  wa' 
heights  and  phases  in  the  numerical  imxlel.  Computations  were  performed  using  tfie 
extended  fiilly-nonlmear  Boussine.si,  model  code  of  Wei  et  al  (1995),  which  is  capabW 
propagatiiig  waves  m  tl.e  large  water  depths  u.sed  in  this  experiment.  Second  and 

•if  ‘•‘’'•‘DMted  l,ased  on  the  entire  experimental  run,  covering 

about  800  wave  periods  Figure  3  shows  the  evolution  of  significant  wave  height  up  tbf 
beach  slope  and  through  the  surf.one,  while  Figure  4  shows  the  evolution  of  skewn^ 

and  asymmetry.  Reproduction  of  measured  values  by  the  numerical  model  is  good 
ootli  cases.  ° 
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Figure  1:  Bottom  configuration  and  wave  gage  locations  for  experiments  of  Masc  and 
LKirby  (1992). 


Figure  2:  Sample  of  measured  (solid)  and  predicted  (dashed)  time  series  of  elevation 
T)  for  Run  2  of  M.isc  and  Kirby  (1992).  Mejisurements  at  h  =  10cm,  corresponding  to 
analysis  of  <lis.sip.ation  rates  below. 


Analysis  of  Data 


Having  vcrific<l  thal  the  nnmerical  model  is  capable  “ . 

Kcond  and  third  moment  statistics  in  the  evolntion  of  a  shoalmg  and  breaking  ran^ 
wave  train,  we  may  now  examine  the  results  tor  energy  d^ipation  m  the  model  ba^ 
on  a  direct  analysis  of  the  laboratory  data.  Figure  6  shows  a  short  s^ment  of  th 
“cord  of  surface  displacement  ,(l)  and  dissipation  «a(t)  at  h  =  10cm.  Di^ipation  ra(0 
L  rrnpnmd  directly  from  the  measured  data  using  (25).  We  compute  the  .»mo<>th^ 
power  spectrum  of  each  of  these  quantities  for  the  entire  data  run,  according  to  the 

definitions 


s,(")  = 
s..{")  = 


2A/ 

2A/ 


(26) 

(27) 


where  ii.  is  the  Fourier  transform  of  the  dissipation  term  and  where  brackets  here 
rnd"atrense!L  averaging.  Results  for  the  Run  2  data  at  h  =  10cm  are  shown  in 

Figure  C. 


f: 


Figure  5i  Time  htetory  of  surface  elevation  and  computed  loss  «a(t)  for  9  seconds  of 
Run  2,  showing  two  strong  breaking  events  at  a  10cm  dept  . 

Returning  to  the  schematic  frequency  domain  model,  we  may  rearrange  the  original 

model  c<iuation  /2g) 

>ln.x  H - -  -QnAn  ' 


into  tlie  form  of  an  energy  equation 


(29) 


m 


0.5  1 


“  filf^) 


Figure  G:  Power  spectrum  S,,{f)  of  surface  displacement  (lower  trace)  and  St{f)  of 
energy  loss  (upper  trace)  for  Mase  and  Kirby  (1992)  Run  2  data. 


The  (juantity  in  brackets  on  the  left  hand  sicle  of  (29)  represents  the  contribution  to 
the  wave  energy  flux  from  the  n’th  frequency  component,  or  F„.  The  quantity  on  the:  V 
right  is  the  contribution  to  the  loss  of  wave  energy  at  that  frequency,  or  C6t»-  Using 
the  definitions  of  power  spectral  densities  (26)  and  (27),  we  may  write  the  dissipation 
coefficient  in  the  form 

„  -  ^  (30) 

Results  similar  to  those  presented  in  Figure  G  may  now  be  \itilizcd  to  determine  the^ 
form  of  o„.  Analysis  of  a  number  of  data  sets  has  indicated  that  the  spectral  tail 
of  Sp(/)  tends  to  have  an  /"^  dependence  on  frequency  after  breaking  is  established 
This  is  consistent  with  the  notion  that  the  wave  form  tends  towards  a  sawtooth  shape, 
with  a  vertical  front  face  and  a  linear  back  slope.  We  note  that  this  is  somewhat 
simplistic  representation  of  the  waves,  as  it  implies  that  the  wavefield  would  have 
significant  asymmetry  and  zero  skewness,  whereas  the  measuretl  data  exhibits  a  bal 
ance  between  skewness  and  asymmetry  in  the  inner  stjrfzone.  Nevertheless,  the  / 
dependence  seems  to  characterize  several  of  the  data  sets  which  have  been  examined 

extremely  well. 

In  contrast.  Figure  6  shows  that  the  dependence  of  5,^(/)  on  /  is  relatively  weak,  so  . 
that  this  quantity  can  be  taken  to  be  constant.  This  is  consistent  with  the  notion  that  ^ 
the  dissipation  term  luts  the  character  of  a  st^iuence  of  Isolated,  spike-like  processes. 
Taken  together,  these  results  indicate  that  the  dominant  frequency  dependence  of  ttnW 

is  given  by  ^ 

or,,  oc  S,,(ti)"'  (31) 


I'  and  wc  may  further  infer  that  this  dependence  translates  into  a  dependence  on  in 
the  inner  surfzone.  Figure  7  shows  a  sample  distribution  of  a{f)  as  calculated  from 
(30),  together  with  a  best  quadratic  fit  to  the  calculated  distribution.  The  fit  is  cen- 
^  tered  fairly  close  to  the  frequency  of  the  spectral  peak  at  IHz,  and  the  dependence  of 
a  on  /  is  clearly  quadratic  towards  higher  frequencies. 


[Figure  7:  Sample  a{f)  deduced  from  data  (according  to  (30))  together  with  a  best 
^quadratic  fit. 

: Effect  of  Choice  F  =  1  on  Modelled  Waves 


The  results  of  the  previous  section  strongly  imply  that  a  value  of  F  =  0  should  be 
chosen  in  the  model  (3)*(5).  CGE  have  shown  that  the  choice  of  F  does  not  introduce 
a  strong  bias  in  the  prediction  of  evolving  power  spectra.  It  appears  that  there  is  a 
preferred  spectral  shape  that  is  obtained  in  shallow  water,  for  which  the  reasons  are 
still  unclear.  The  effect  of  redistributing  the  loss  differently  across  the  spectrum  serves 
mainly  to  enhance  or  suppress  the  nonlinear  transfer  of  energy  needed  to  maintain 
the  target  spectral  shape. 

In  particular,  the  choice  F  =  1  implies  that  the  rate  of  energy  loss  is  the  same  in 
>all  spectral  components,  which  accounts  for  the  entire  pattern  of  overall  energy  loss 
in  evolving  waves  (Mase  and  Kirby,  1992;  Eldeberky  and  Battjes,  199G).  However,  if 
'  the  dissipation  term  is  chosen  to  account  for  all  changes  in  spectral  energy  density, 
there  is  necessarily  a  parallel  suppression  of  all  nonlinear  energy  transfer  across  the 
spectrum  during  the  breaking  process.  This  loss  of  an  active  transfer  should  suppress 
the  imaginary  part  of  the  bispectrum,  and  would  be  evidenced  by  a  loss  of  front-to- 
back  asymmetry  in  the  wave  form,  or  a  loss  of  statistical  asymmetry.  CGE  have  shown 
the  consequence  of  choosing  F  =  1  on  the  prediction  of  third  moment  statistics,  which 


;irc  scviM.  Iv  A  iihhc  .•x.iiiipic  i.s  |)H*vi<lo«l  liy  Liu <(11)90,  see  pages  55^ 

v7),  wlio  modelled  tlie  Lieakiiig  decay  ol  a  periodic  wave  using  a  dissipation  coefficient 
wlii(  li  w<Ls  uiiilorni  across  tli<“  imxlelled  sp(!ctruni.  'I'lie  ri^sulting  wave  crests  in 
snrl/one  clearly  lose  their  asymmetry,  which  provi<les  direct  evidence  that  nonlincS 
energy  transler  across  the  spectrum  has  been  siijipressed. 

In  contrjust,  the  choice  /•’  =  0  strongly  concentrates  dissipation  at  higher  freqtie^ 
cies.  As  a  residt,  there  is  necessarily  a  high  rate  of  noidinear  transfer  of  energy  ac 
tin;  sfiectrum  licmi  low  to  high  fre<|iieiicies,  in  oder  to  maintain  the  tail  of  the  sp 
trum.  This  accounts  for  the  presence  of  a  strong,  negative  asyininetry.  This  choice 
iniich  more  sensible  in  light  of  our  understanding  of  the  nonlinear  processes  going  oS 
in  this  legion,  cuid  agrees  with  a  direct  analysis  of  the  behavior  of  a  well-tested  tiin* 
domain  model  of  the  breaking  process. 

Conclusions 

1  he  I  (‘.suits  of  this  study  indicat(*  that  di.ssi|>ation  due  to  wave  breaking  shoiUd 
be  biiiscd  strongly  towards  higher  frtHjuencies  in  spectral  calculations,  and  thati^M 
f  dependente  in  the  dissipation  coeflicent  come.s  closest  to  matching  the  desiri^ 
structure  of  the  breaking  terms  using  a  .simple  functional  dependence  This  resul' 
supports  the  conclu.sions  of  Chen  et  al  (1990)  and  is  obtained  by  an  entirely  differen 
rout(‘,  Ixung  based  on  an  aidysis  of  l.h(?  f r(‘(|U('ncy  .structure!  ol  tlu!  dis.sipation  term  in 
a  tiiiK'-domain  wave  (‘volution  model.  «! 

A  mu(  h  mor(‘  compr(‘h(.‘n.siv(‘  d(‘.scription  of  th(‘  r(!.sidt.s  of  this  study  i.s  currently 
being  pi(!p,ii(‘d,  and  will  consid(!r  th(‘  form  ol  tli(‘  dis.sipation  coc'flicient  in  the  outiu', 
snrf/one  as  w(‘ll  as  m  th(‘  (‘.stablish(‘d  inn(‘r  surf/oue.  In  i)articular,  there  is  so'mn; 
indication  that,  in  rc‘gions  where  the  breaking  events  are  infrequent,  that  there  is^i' 
tendency  for  o„  to  b(‘  slightly  negativ(!  at  fr(‘(|uenci(‘.s  b(‘low  the  wind  wave  peak.  TiiS 
result  coidd  b(‘  thought  of  as  th(‘  effect  of  distributing  a  s(‘t  of  localized  momentniA 
sourc(.‘.s  in  th(‘  domain.  L,ach  breaking  wav(‘  (‘V(!nt  woidd  impart  a  kick  to  the  waiter 
column  in  the  mann(‘r  (le.scrib(‘d  by  Happ  and  Melville  (1990).  The  irregular,' widef 
spaced-in-time  nature  of  these  kicks  translates  into  a  wave  generating  mcchanism^i^ 
low  frc(|uency,  and  could  contribute  to  the  growth  of  the  low-frequency  wave  clima^ 
in  the  surfV.on(‘.  , 
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